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We present a search for the decay B + —* t + v using 383 x 10 6 BB pairs collected at the T(45) 
resonance with the BABAR detector at the SLAC PEP-II B Factory. We select a sample of events 
with one completely reconstructed tag B in a hadronic decay mode (B~ — > D^*'° X~), and examine 
the rest of the event to search for a B + — > t + v decay. We identify the r lepton in the following 
modes: r + — * e + iSU, r + — » p^vV, r + — > tt + V and r + — + 7r + 7r°77. We find a 2.2 a excess in 
data and measure a branching fraction of B(B + -> t + v) = (1.8t^(stat.) ± 0.4(bkg. syst.) ± 
0.2(other syst.)) x 10 -4 . We calculate the product of the B meson decay constant /b and \V u b\ to 
be f B ■ \V ub \ = (lO.lt^(stat.)^(syst.)) x HT 4 GeV 
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PACS numbers: 13.20.-v, 13.25.Hw 



The study of the purely leptonic decay B + — > t+v [l[ is 
of particular interest because it is sensitive to the prod- 
uct of the B meson decay constant f B , and the abso- 
lute value of Cabibbo-Kobayashi-Maskawa matrix ele- 
ment V ub HQ. In the Standard Model (SM), the decay 
proceeds via quark annihilation into a W+ boson, with 
a branching fraction given by: 



B{B J[ 



8?r 



1 



Tb+Ib\ V ui 



(1) 

where Gf is the Fermi constant, t b + is the B + life- 
time, and m B and m T are the B + meson and r lep- 
ton masses. Using \V ub \ = (4.31 ± 0.30) x 10~ 3 from 
experimental measurements of semileptonic B decays Q 
and f B = 0.216 ± 0.022 GeV from lattice QCD @, the 
SM estimate of the branching fraction for B + — > t+v is 
(1.5 ±0.4) x 10" 4 . 

The process B + — > t+v is also sensitive to extensions 
of the SM. For instance, in two-Higgs doublet models [6j 
and in the MSSM [7|, l8J it could be mediated by charged 
Higgs bosons. The branching fraction measurement can 
therefore also be used to constrain the parameter space 
of extensions to the SM. 

The B + — * /i+v and B + — * e+v decays are signif- 
icantly helicity suppressed with respect to the B + — > 
t+v channel. However, a search for B + — > t+v is 
experimentally more challenging, due to the presence 
of multiple neutrinos in the final state, which makes 
the experimental signature less distinctive. In a pre- 
viously published analysis using a sample of 383 x 10 6 
T(AS) — > BB decays, based on the reconstruction of a 
semileptonic B decay on the tag side, the BaBAR col- 
laboration set an upper limit B(B + — > t+v) < 1.7 x 
10~ 4 at the 90% confidence level (CL) [§]. The Belle 
Collaboration has reported evidence from a search for 
this decay and the branching fraction was measured to 
be B(B+ -> t+v) = (1.79i°-^(stat.)io5?(syst.)) x 10~ 4 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II storage ring. The sample 
corresponds to an integrated luminosity of 346 fb _1 at the 
T(AS) resonance (on-resonance) and 36.3 fb _1 taken at 
40MeV below the T(4S) resonance (off-resonance). The 
on-resonance sample contains 383 x 10 6 BB decays. The 
detector is described in detail elsewhere [ll[. Charged- 
particle trajectories are measured in the tracking sys- 
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tern composed of a five-layer silicon vertex detector and 
a 40-layer drift chamber (DCH), operating in a 1.5 T 
solenoidal magnetic field. A Cherenkov detector is used 
for 7T— K discrimination, a Csl calorimeter for photon de- 
tection and electron identification, and the flux return of 
the solenoid, which consists of layers of steel interspersed 
with resistive plate chambers or limited streamer tubes, 
for muon and neutral hadron identification. 

In order to estimate signal selection efficiencies and to 
study physics backgrounds, we use a BABAR Monte Carlo 
(MC) simulation based on GEANT4 [12]. In MC simu- 
lated signal events one B + meson decays to t+v and the 
other into any final state. The BB and continuum MC 
samples are, respectively, equivalent to approximatively 
three times and 1.5 times the accumulated data sample. 
Beam-related background and detector noise are taken 
from data and overlaid on the simulated events. 

We reconstruct an exclusive decay of one of the B 
mesons in the event (tag B) and examine the remaining 
particle(s) for the experimental signature of B + — > t+v. 
In order to avoid experimenter bias, the signal region 
in data is blinded until the final yield extraction is per- 
formed. 

The tag B candidate is reconstructed in the set of 
hadronic B decay modes B~ — > D^^X^ jlj, where 
X~ denotes a system of charged and neutral hadrons 
with total charge —1 composed of n\Tt+- , nzK+ , n^K^, 
n^ir , where n\ + ri2 < 5, < 2, and 714 < 2. We re- 
construct D*° -> D°n ,D Q "/; D° -► K~tt + , R-tt+tt , 
K~ir + ir~ir + , K®ir + Tr~ and Kg — > ir + ir~ . The kinematic 
consistency of tag B candidates is checked with the beam 
energy-substituted mass rngs = — p B and the en- 

ergy difference AE = E B — y/s/2. Here y/s is the total 
energy in the T(AS) center-of-mass (CM) frame, and p B 
and E B denote, respectively, the momentum and energy 
of the tag B candidate in the CM frame. The resolution 
on AE is measured to be ctab = 10 — 35 MeV, depending 
on the decay mode; we require \AE\ < 3a ae- The purity 
V of each reconstructed B decay mode is estimated, using 
on-resonance data, as the ratio of the number of peaking 
events with m^s > 5.27 GeV/ c 2 to the total number of 
events in the same range. If multiple tag B candidates 
are reconstructed, the one with the highest purity V is 
selected. If more than one candidate with the same pu- 
rity is reconstructed, the one with the lowest value of 
\AE\ is selected. From the dataset obtained as described 
above, we consider only those events in which the tag B 
is reconstructed in the decay modes of highest purity V . 
The set of decay modes used is defined by the require- 
ment that the purity of the resulting sample is not less 
than 30%. 

The background consists of e + e~ —^qq(q — u, d, s, c) 
events and other T(AS) — * B°B° or B + B~ decays in 
which the tag B candidate is mis-reconstructed using par- 
ticles coming from both B mesons in the event. To reduce 
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the e + e _ — > background, we require | cos6>^ B | < 0.9, 
where 6^, B is the angle in the CM frame between the 
thrust axis [HI of the tag B candidate and the thrust 
axis of the remaining reconstructed charged and neutral 
candidates. 

In order to determine the number of correctly recon- 
structed B + decays, we classify the background events 
in four categories: e + e~ — > cc; e+e~ — > uu,dd,ss; T(AS) 
-> B°B°; and T(45) -> B+B~ . The to es shapes of 
these background distributions are taken from MC simu- 
lation. The normalization of the e + e~ — > cc and e + e~ — > 
mm, dd, ss backgrounds is taken from off-resonance data, 
scaled by the luminosity and corrected for the different 
selection efficiencies evaluated with the MC. The normal- 
ization of the B°B Q , B + B~ components are obtained by 
means of a x 2 fit to the ibes distribution in the data side- 
band region (5.22GeV/c 2 < m ES < 5.26GeV/c 2 ). The 
number of background events in the signal region (toes > 
5.27GeV/c 2 ) is extrapolated from the fit and subtracted 
from the data. We estimate the total number of tagged 
B's in the data to be N B = (5.92±0.11(stat)) x 10 5 . Fig- 
ure [1] shows the tag B candidate toes distribution, with 
the combinatorial background, estimated as the sum of 
the four components described above, overlaid. 



x10 3 




m ES (GeV/c 2 ) 

FIG. 1: Distribution of the energy substituted mass, jties, of 
the tag B candidates in data. The combinatorial background 
is overlaid. 



After the reconstruction of the tag B meson, a set of 
selection criteria is applied to the rest of the event (recoil) 
in order to enhance the sensitivity to B + — > t + v decays. 
We require the presence of only one well-reconstructed 
charged track (signal track) with charge opposite to that 
of the tag B. The signal track is required to have at least 
12 hits in the DCH, momentum transverse to the beam 
axis, pTi greater than 0.1 GeV/c, and the point of closest 
approach to the interaction point less than 10 cm along 
the beam axis and less than 1.5 cm transverse to it. 

The r lepton is identified in four decay modes con- 
stituting approximately 71% of the total t decay width: 
t + — * e + vV, t + — > p + vV, t + — > 7r + F, and r + — > 7r + 7r°77. 
Particle identification criteria on the signal track are used 
to separate the four categories. The r + — > 7r + 7r°77 sample 



is obtained by associating the signal track, identified as 
pion, with a 7r° reconstructed from a pair of neutral clus- 
ters with invariant mass between 0.115 and 0.155 GeV/c 2 
and total energy greater than 250 MeV. In case of multi- 
ple 7T + 7r° candidates, the one with largest center-of-mass 
momentum p* . is chosen. 

We place a mode-dependent cut on | cos 9^ B | to reduce 
the background due to continuum events and incorrectly 
reconstructed tag B candidates (combinatorial). The re- 
maining sources of background consists of B + B~ events 
in which the tag B meson was correctly reconstructed 
and the recoil contains one track and additional parti- 
cles that are not reconstructed by the tracking detectors 
and calorimeter. MC simulation shows that most of this 
background is from semileptonic B decays. 

We define the discriminating variable -E cx tra as the sum 
of the energies of the neutral clusters not associated with 
the tag B or with the signal ir° from the r + — > tt + tt V 
mode, and passing a minimum energy requirement. The 
required energy depends on the selected signal mode and 
on the calorimeter region involved and varies from 50 
to 70 MeV. Signal events tend to peak at low -E extra 
values, whereas background events, which contain addi- 
tional sources of neutral clusters, are distributed toward 
higher E cxtr a values. 

Other variables used to discriminate between signal 
and background are the CM momentum of the signal 
candidates, the multiplicities of low px charged tracks 
and of 7T° candidates in the recoil, and the direction of 
the missing momentum four-vector in the CM frame. For 
the t + — > 7r + 7r°77 mode, we exploit the presence of the 7r° 
in the final state and the dominance of the decay through 
the p + resonance by means of the combined quantity 
x p = [(TO 7r+7r o-TOp)/(r p )] 2 + [(TO 77 -TO w o)/(CT 7r o)] 2 , where 
m^+^o is the reconstructed invariant mass of the 7r + 7r° 
candidate, the reconstructed invariant mass of the 

7T° candidate, m p and V p are the nominal values [H for 
the p mass and width, m T o is the nominal ir° mass and 
cr n o = 8 MeV/ c 2 is the experimental resolution on the ir° 
mass determined from data. 

We optimize the selection by maximizing s/\Js + b us- 
ing the B + B~ MC and signal MC, where b is the ex- 
pected background from B + B~ events and s is the ex- 
pected number of signal events in the hypothesis of a 
branching fraction of 1 x 10" 4 . The optimization is per- 
formed separately for each r decay mode and with all the 
cuts applied simultaneously in order to take into account 
any correlations among the discriminating variables. The 
optimized signal selection cuts are reported in Table |TJ 

We compute the signal selection efficiency as the ratio 
of the number of signal MC events passing the selection 
criteria to the number of signal events that have a cor- 
rectly reconstructed tag B candidate in the signal region 
toes > 5.27GeV/c 2 . We evaluate the efficiencies on a sig- 
nal MC sample which is distinct from the sample used in 
the optimization procedure. A small cross-feed in some 
modes is estimated from MC and is taken into account 
in the computation of the total efficiency. 
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TABLE I: Optimized selection criteria for each r decay mode. 



Variable 



f-i 



E, 



extra 



( GeV) < 0.160 < 0.100 < 0.230 



7T° multiplicity 

Track multiplicity 1 

\cos6* TB \ < 0.9 

Pt * rk (GeV/c) <1.25 

cos8* miaa < 0.9 
pXo(GeV/c) 

X p 

E^o (GeV) - 




1 

< 0.9 
< 1.85 



< 2 

< 2 

< 0.7 
> 1.5 

< 0.5 



< 0.290 



1 

< 0.7 

< 0.55 
> 1.5 

< 2.0 
> 0.250 



TABLE III: Observed number of on-resonance data events 
in the signal region compared with the number of expected 
background events. 

t decay mode Expected background Observed 



r + — > e + vV 


1.47 ± 1.37 


4 


T + — * /l + l/J7 


1.78 ± 0.97 


5 


T + — > 7T + 77 


6.79 ± 2.11 


10 


r + — > 7r + 7r°77 


4.23 ± 1.39 


5 


All modes 


14.27 ± 3.03 


24 



The total efficiency for each selection is given by: 



(2) 



where is the efficiency of the selection i for the r decay 
mode j, ridoc = 7 is the number of r decay modes that 
can contribute to the reconstructed modes and fj are 
the fractions of the r decay mode as estimated from the 
signal MC sample with a reconstructed tag B. Table ITT1 
shows the estimated efficiencies. 



TABLE II: Efficiency (in percent) of the most relevant r de- 
cay modes (rows) to be selected in one of the four modes con- 
sidered in this analysis (column). The All decay row shows 
the selection efficiency of each reconstruction mode, adding 
the contribution from the previous rows, weighted by the de- 
cay abundance at the tag selection level fj. The last row 
shows the total signal selection efficiency. The uncertainties 
are statistical only. 



Mode 



/' 



Other 



19.3 ± 










.1 

10.8 ± 0.9 

0.1 ± 0.1 







0.4 ± 0.2 
1.3 ± 0.3 
19.7 ± 1.3 

1.5 ± 0.2 


0.2 ± 0.1 
0.3 ± 0.2 






0.5 ± 0.2 

7.0 ± 0.5 



1.8 ± 0.4 

0.1 ± 0.1 



All dec. 



3.1±0.2 1.7±0.1 



2.9±0.2 2.2±0.2 



Total: 



9.8 ± 0.3 



To determine the expected number of background 
events in the data, we use the final selected data sam- 
ples with .Eextra between and 2.4 GeV. We first 
perform an extended unbinned maximum likelihood fit 
to the niES distribution in the E ext rn sideband region 
0.4 GeV < Scxtra < 2.4 GeV of the final sample. For the 
peaking component of the background we use a probabil- 
ity density function (PDF) which is a Gaussian function 
joined to an exponential tail (Crystal Ball function) [TJ. 



As a PDF for the non-peaking component, we use a 
phase space motivated threshold function (ARGUS func- 
tion) [15j . From this fit, we determine a peaking yield 
-^pfc de data aim signal shape parameters, to be used in later 
fits. We apply the same procedure to B + B~ MC events 
which pass the final selection and determine the peaking 



yield N\ 



side.MC 
pk 



To determine the MC peaking yield in 
the #extra signal region N^ k 6,MC , we fit to E s in the E cxtla 
signal region of the B + B~ MC sample with the Crys- 
tal Ball parameters fixed to the values determined in the 
-Eextra sideband fits described above. Analogously, we fit 
the rriES distribution of data in the -Eextra signal region 
to extract the combinatorial background n com b, evalu- 
ated as the integral of the ARGUS shaped component 
in the rriES > 5.27 GeV/ c 2 region. We estimate the total 
expected background in the signal region as: 



b = 



N' 



sig,MC 



pk 



sido,MC 



x Nl 



side, data 
pk 



^comb- 



(3) 



pk 



After finalizing the signal selection criteria, we measure 
the yield of events in each decay mode in on-resonance 
data. Table Hill reports the number of observed events to- 
gether with the expected number of background events, 
for each r decay mode. Figure O shows the E extYa distri- 
bution for data and expected background at the end of 
the selection. The signal MC, normalized to a branching 
fraction of 3 x 1CU 3 for illustrative purposes, is overlaid 
for comparison. The E ey ± ra distribution is also plotted 
separately for each r decay mode. 

We combine the results on the observed number of 
events rii and on the expected background hi from each 
of the four signal decay modes {n c h) using the esti- 
mator Q = C(s + b)/C(b), where C(s + b) and C{b) are 
the likelihood functions for signal plus background and 
background-only hypotheses, respectively: 



i=l 



i=l 



(4) 

The estimated number of signal candidates Si in data, for 
each decay mode, is related to the B + — > t + v branching 
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e Data 

| total background prediction 
1 combinatorial background 
Signal MC (scaled to BF=3x10 3 ) 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 
E e x,ra(GeV) 




0.5 1 1.5 2 



- «l) 









E extra (GeV) 



FIG. 2: -Eoxtra distribution after all selection criteria have 
been applied. The upper plot shows the distribution of all the 
modes combined while lower plots show the (a) r + — » e + vv, 
(b) t + — > /i + ^77, (c) r + — » 7r + 77, and (d) r + — + 7r + 7r°I' modes 
separately. The on-resonance data (black dots) distribution is 
compared with the total background prediction (continuous 
histogram). The hatched histrogram represents the combina- 
torial background component. B + — * t + v signal MC (dashed 
histogram), normalized to a branching fraction of 3 x 10 -3 for 
illustrative purposes, is shown for comparison. 



fraction by: 



_tag J "B+ 



(5) 



where N^f is the number of tag B + mesons correctly 
reconstructed, e^ s and e s ?f are the tag B efficiencies in 
generic BB and signal events respectively, and are the 
signal efficiencies defined in equation [2] We fix the ratio 
e^/eg & = 0.939 ± 0.007(stat.) to the value obtained 
from MC simulation. 

We estimate the branching fraction (including statisti- 
cal uncertainty and uncertainty from the background) by 
scanning over signal branching fraction hypotheses and 
computing the value of C{s + b)/C(b) for each hypothe- 
sis. The branching fraction is the hypothesis which mini- 
mizes the likelihood ratio— 21nQ = — 2\n(£(s + b)/ C(b)), 
and we determine the statistical uncertainty by finding 
the points on the likelihood scan that occur at one unit 
above the minimum. 

The dominant uncertainty on the background predic- 



tions bi is due to the finite B + B~ MC statistics. We 
also check possible systematic effects in the estimation 
of combinatorial background by means of a sample of 
events with looser selection requirements; we find it to 
be negligible with respect to the statistical uncertainty. 
The background uncertainty is incorporated in the like- 
lihood definition used to extract the branching fraction, 
by convolving it with a Gaussian function with standard 
deviation equal to the error on bi [16l |. 

The other sources of systematic uncertainty in the de- 
termination of the B + — > t + v branching fraction come 
from the estimation of the tag yield and efficiency and 
the reconstruction efficiency of the signal modes. We es- 
timate the systematic uncertainty on the tag B yield and 
reconstruction efficiency by varying the MC B + B~ non- 
peaking component of the toes shape, assigning a sys- 
tematic uncertainty of 3% on the branching fraction. The 
systematic uncertainties due to mismodeling of charged 
particle tracking efficiency, E ext ra. shape, particle identifi- 
cation efficiency, ir° reconstruction and signal MC statis- 
tics depend on the r decay mode. The uncertainty on the 
branching fraction is evaluated for each mode separately. 
We obtain the total contributions due to tracking and 
-E'extra systematics by adding linearly the contributions 
of each decay channel. The total contributions due to 
MC statistics and particle identification are obtained by 
adding systematics uncertainties of each reconstruction 
mode in quadrature. 

We check the low pt charged track multiplicity dis- 
tribution agreement between data and MC with a sam- 
ple enriched in background by loosening the selection 
criteria. The disagreement, which is mode dependent, 
is quantified by comparing the MC PDF with the data 
PDF. We correct the MC to reproduce the distribution 
in data and apply the correction to the signal MC distri- 
bution. We take 100% of the correction as a systematic 
uncertainty, resulting in a total systematic uncertainty of 
5.8% on the branching fraction. 

The systematic uncertainty due to the E cyAva mismod- 
eling is determined by means of a data sample containing 
events with two non-overlapping tag B candidates. The 
sample is selected by reconstructing a second B meson 
in a hadronic decay mode B~ — > flW X~ on the recoil 
of the tag B. In addition to the requirements on the 
tag B described above, we consider only second B candi- 
dates satisfying \ AE\ < 50MeV and m ES > 5.27GeV/c 2 
having opposite charge to that of the tag B. If multiple 
candidates are reconstructed, the one with the highest 
purity V is selected. We compare the distribution of the 
total energy of the unassigned neutral clusters -Eextra in 
data and in MC. We compute the ratio of the number of 
events in the signal region of each r mode to the total 
number of events in the sample. For each r mode, we 
evaluate the systematic uncertainty, comparing the ra- 
tio estimated from MC to the ratio estimated from data. 
This procedure results in a 8.8% systematic uncertainty 
on the branching fraction. Table IIVI shows the contribu- 
tions in percent to the systematic uncertainties on the 
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TABLE IV: Contributions (in percent) to the systematic un- 
certainty on the branching fraction due to signal selection 
efficiency for different selection modes. 



Source of systematics 


e+ 


M + 


7T+ 




Total 


MC statistics 


3.1 


0.6 


1.5 


2.6 


4.3 


Particle Identification 


1.5 


1.3 


0.2 


0.2 


2.0 










1.4 


1.4 


Tracking 


3.7 


0.4 


0.1 


1.6 


5.8 


-^extra 


4.7 


0.6 


0.9 


2.6 


8.8 


Signal B 










11.6 


Tag B 










3 


Total 










12 



branching fraction. 

In summary, we measure the branching fraction 

B(B+ -> t+v) = (1.8±g;| ± °- 4 ± 0-2) x (6) 

where the first error is statistical, the second is due to the 
background uncertainty, and the third is due to other sys- 
tematic sources. Taking into account the uncertainty on 
the expected background, as described above, we obtain 
a significance of 2.2 a. 

Using Eq. [TJ we calculate the product of 



the B meson decay constant fg and \V u b\ to be 
!b ■ \Vub\ = (10.ll|j(stat.)t^(syst.)) x lfr 4 GeV. We 
also measure the 90% C.L. upper limit using the CL S 
method [H to be B(B+ t+v) < 3.4 x lfr 4 . 

The combination of this measurement with the BABAR 
result obtained using semileptonic tags, based on a sta- 
tistically independent data sample, and reported in 0], 
yields: 

B(B+ - t+v) = (1.2±0.4 stat ±0.3 bkg ±0.2 syst .)xl0- 4 . 

The significance of the combined result is 2.6 a including 
the uncertainty on the expected background (3.2 a if this 
uncertainty is not included). 
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